A 


RADC-TR-89-217 
Rnal  Technical  Report 
Octobar  1989 


IONOSPHERIC  HEATING  BY 
OBLIQUELY  TRANSMITTED  HIGH 
POWER  RADIO  WAVES 


Mas;$achusetts  Institute  of  Technology 


Ming*Chang  Lee 


APPROVED  FOR  PUBUC  RELEASE;  DISTRIBUTION  UNLIMITED. 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  Systems  Command 
Griffiss  Air  Force  Base,  NY  13441-5700 


?  0  y  ^ 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Division  (PA) 
and  Is  releasable  to  the  National  Technical  Infonnation  Services  (NTZS)  At 
NTIS  It  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR~89-217  has  been  reviewed  and  is  approved  for  publication. 


APPROVED:  (P 

STANFORD  P.  YUKON 
Project  Engineer 


APPROVED: 

JOHN  K.  SCHINDLER 
Director  of  Electromagnetics 


FOR  THE  COMMANDER: 

JOHN  A.  RITZ 

Directorate  of  Plans  &  Programs 


If  your  address  has  changed  or  If  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  If  the  addressee  Is  no  longer  employed  by  your 
organization,  please  notify  RADC  (  EECP  )  Griffiss  AFB  NY  13441-5700. 
This  will  assist  us  in  maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or 
notices  on  a  specific  document  require  that  It  be  returned. 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
0MB  No.  0704-0188 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

N/A 

2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

N/A 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

N/A 

6a.  NAME  OF  PERFORMING  ORGANIZATION 
Massachusetts  Institute 
of  Technology 

6b.  OFFICE  SYMBOL 
at  tppIkabM) 

6c  ADDRESS  (Oty,  State,  and  ZIP  Coda) 

Plasma  Fhslon  Center 

Cambridge  MA  02139 

8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Rome  Air  Development  Center 

Bb.  OFFICE  SYMBOL 

Of  appBttbM) 

EECP 

8c  ADDRESS  (Oty,  State,  j/id  ZIP  Code) 

Hanscom  AFB  MA  01731-5000 

1b.  RESTRICTIVE  MARKINGS 

N/A 


3.  DISTRIBUTION /AVAILABIUTY  OF  REPORT 
Approved  for  public  release; 
distribution  unlimited. 


$.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

RADC-TR-89-217 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Rome  Air  Development  Center  (EECP) 


7b.  ADDRESS  (Oty,  Statw,  and  Z/P  Coda) 


Hanscom  AFB  MA  01731-5000 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F30602-81-C-0193 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

62702F 


PROJECT 

TASK 

NO 

NO 

4600 

16 

WORK  UNIT 
ACCESSION  NO 


11.  TITLE  (Indum  Saepriiy  Oaadfeadon) 


IONOSPHERIC  HEATING  BY  OBLIQUELY  TRANSHirrED  HIHG  POWER  RADIO  WAVES 


1Z.  PERSONAL  AUTHOR(S) 

Mln-Chang  Lee 


13a.  TYPE  OF  REPORT  113b.  TIME  COVERED  114.  DATE  OF  REPORT  (Yaar.  IMomfi,  Day)  1 15.  PACE  COUNT 

Final  I  from  Feb  88  to  Dec  88  f  October  1989  I  36 


16.  SUPPLEMENTARY  NOTATION 

N/A 


IT.  ■  •  COSATI  COOES  IB.  SUBJECT  TERMS  (Cor»f<ii«a  on  ravaraa  if  nacasary  and  kMnofy  by  btodt  numbtri 

FIELD  GROUP  SUB-GROUP 

20  14  Oblique  Ionospheric  Modification 


19.  ABSTRACT  (ConcfniM  on  raaoraa  if  ntnttarf  and  idtmif^f  by  bfotk  nombar) 

Anisotropic  effects  of  Ionospheric  Irregularities  on  Faraday  polarization  fluctuations 
have  been  investigated.  It  was  found  that  sheet  Ilka  Irregularities  are  more  effective 
than  rod  like  Irregularities  In  causing  Faraday  polarization  fluctuations.  Measurement 
of  Che  fluctuation  spectrum  may  be  used  Co  diagnose  Ionospheric  Irregularities.  Non¬ 
linear  scattering  of  HF  waves  by  ionospheric  irregularities  has  also  been  investigated. 

It  was  shown  that  HF  ionospheric  heating  can  be  further  enhanced  by  nonlinear  scattering 
of  VLF  waves.  Artificial  spread  F  caused  by  ionospheric  heating  with  either  vertically 
or  obliquely  incident  high  power  HF  waves  has  been  investigated  by  ray-tracing  simulation 
’  codes. 


20.  DISTRIBUTION /AVAILABIUTY  OF  ABSTRACT 
□  UNCLASSJFICORJNUMITEO  □  SAME  AS  RPT. 


22«.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Stanford  P.  Yukon 


DO  Form  1473,  JUN  86 


121.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22b.  TELEPHONE  (Jndwdb  ArMCodb)  22c  OFFICE  SYMBOL 

(617)  377-2985 


PlwJouStdrnomarRObloMM.  SECURITY  CLASSIFICATION  OF  THIS 

UNCLASSIFIED 


Final  Report  to  Rome  Air  Development  Center  for  Sponsored  Research  in  “Ionospheric 
Heating  by  Obliquely  Transmitted  High  Power  Radio  Waves  ”  during  February  3  - 
December  3,  1988.  Prepared  by 


Dr.  Min-Chang  Lee 
Principal  Investigator 


1.  Introduction 

Continued  research  under  the  sponsorship  of  RADC  hcus  been  actively  conducted  at 
Massachusetts  Institute  of  Technology  in  ionospheric  heating  by  obliquely  transmitted  high 
power  radio  waves.  During  the  period  of  February  3  -  December  3,  1988,  research  effort 
centered  on  (1)  analysis  of  induced  ionospheric  effects  ,  (2)  diagnosis  of  oblique  heating 
effects  by  ionosondes,  and  (3)  planning  of  the  ECRS-Bermuda  oblique  heating  experiment 
with  RADC. 

Because  of  the  oblique  wave  propagation  in  the  ionosphere,  resonant  wave-plasma 
interactions  which  require  strict  frequency  and  wave  vector  matching  relations  are 
not  expected  to  occur.  These  resonant  wave-plasma  interaction  processes,  which 
can  generate  short-scale  (typically,  meter-scale)  ionospheric  density  irregularities  play 
important  roles  in  the  vertical  heating  experiments.  It  is,  however,  expected  that  large- 
scale  (typically,  hundred  meter-scale)  ionospheric  density  irregularities  are  able  to  be 
excited  by  thermal  filzunentation/ self-focusing  instabilities.  These  large-scale  instabilities 
can  produce  prominent  ionospheric  effects  along  the  wave  path  in  a  large  ionospheric 
region.  In  the  following,  we  will  briefly  report  our  analysis  of  these  heater  wave-induced 
ionospheric  disturbances  and  their  effects  on  radio  propagation. 

2.  Heater  wave-induced  ionospheric  irregularities 

Ionospheric  density  irregularities  generated  by  filamentation  instabilities  can  have  two 
independent  wave- vector  directions.  The  irregularities  excited  by  the  o  mode  heater  waves 
are  field  aligned  and  have  wave-vectors  pointing  to  the  direction  perpendicular  to  the 
meridian  plane.  By  contrast,  thf  irregularities  excited  by  the  x  mode  heater  waves  are 
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polarized  in  the  meridian  plane  and  are,  in  general,  not  field  aligned.  The  field  aligned 
nature  of  induced  ionospheric  irregularities,  however,  is  eventually  established  to  reduce 
t  he  diffusion  damping  along  the  geomagnetic  field. 

.4s  illustrated  in  Figure  1,  the  radio  signals  transmitted  vertically  from  an  ionosonde 
can  have  two  different  types  of  ray  trajectories  in  the  presence  of  ionospheric  density 
irregularities.  The  reflection  heights  of  these  ray  trajectories  are  significantly  different  from 
those  of  ray  trajectories  in  the  absence  of  heater  induced  ionospheric  irregularities.  The 
fluctuation  of  reflection  heights,  then,  gives  rise  to  the  spread  F  echoes  on  the  ionograms. 
It  can  be  shown  that  irregularities  polarized  within  the  meridian  plane  can  cause  spread 
F,  while  those  polarized  in  planes  perpendicular  to  the  meridian  planes  can  not. 

We  have  developed  ray  tracing  codes  to  simulate  the  artificial  spread  F  caused  by 
ionospheric  heating  w-ith  either  vertically  or  obliquely  transmitted  high  power  radio  waves. 
.4  linear  ionospheric  density  profile  is  adopted  in  the  ray  tracing  codes.  Further,  we 
also  develop  a  ray  tracing  code  for  radial  density  profile  which  represents  the  typical 
inhomogeneous  laboratory  plasmas.  This  work  is  part  of  our  planned  laboratory  study 
of  nonlinear  wave-plasma  interactions  which  can  simulate  the  oblique  ionospheric  heating 
experiments.  Two  students,  Dan  Moriarty  and  Larry  McKay,  are  working  with  Dr.  Lee 
on  ray  tracing  study  of  radio  propagation  in  perturbed  ionosphere  for  their  thesis  research. 
The  results  of  the  research  will  be  reported  later. 

3.  Anisotnipic  effects  of  ionospheric  irregularities 

The  anisotropic  effects  of  ionospheric  irregularities  on  Faraday  polarization  fluctuations 
have  been  investigated  by  Dr.  Lee  with  a  student,  S.V.  Nghiem.  This  work  extends 
the  earlier  theory  of  Lee  et  al  1982j,  explaining  the  Faraday  polarization  fluctuations  of 
satellite  beacon  signals  at  136  MHz.  In  the  previous  work,  isotropic  correlati'  i  functions 
were  assumed  for  the  ionospheric,  density  fluctuations  for  simplicity.  However,  ionospheric 
density  irregularities  have  prominent  anisotropic  nature,  namely,  they  '^.r.ngate  along  the 
earth's  magnetic  field  to  reduce  diffusion  damping.  Rod-like  and  sKeet-like  cire  the  two 
forms  of  ionospheric  irregularities  commonly  seen  in  the  ionosphere.  Sheet-like  ionospheric 
irregularities  are  observed  in  the  auroral  oval,  while  rod-like  irregularities  at  higher  and 
lower  latitudes.  The  different  geometric  effects  of  these  ionospheric  irregularities  can  be 
exhibited  in  the  scintillation  measurements  of  radio  signals. 

Striking  anisotropic  effects  of  ionospheric  irregularities  are  found  on  the  Faraday 
polarization  fluctuations  of  linearly  polarized  radio  waves,  propagating  along  the  earth's 
magnetic  field.  .Nevertheless,  if  the  wave  propagation  angle  is  not  small  (say,  greater 
than  •")  degrees)  the  anisotropic  effects  of  r^  d-like  irregularities  weaken  significantly  while 
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those  of  sheet-like  irregularities  remain  prominent.  In  general,  the  sheet-like  ionospheric 
irregularities  are  more  effective  than  the  rod-like  irregularities  in  causing  the  Faraday 
polarization  fluctuation  of  radio  waves. 

The  present  work  can  be  baised  on  to  develop  a  radio  technique  for  diagnosing  the 
ionospheric  irregularities  induced  by  the  obliquely  transmitted  radio  waves.  As  shown 
in  Figure  2,  a  linearly  polarized  diagnostic  radio  wave  with  a  slightly  less  frequency  is 
transmitted  along  the  same  path  of  the  heater  wave.  Then,  the  heater  wave-induced 
ionospheric  irregularities  will  cause  the  polarization  fluctuations  of  the  diagnostic  wave. 
Further,  the  nature  of  the  induced  ionospheric  irregularities  can  also  be  deduced  on  the 
basis  of  the  present  theory.  A  preprint  of  the  paper  is  attached  as  .Appendix  .A  of  this 
report . 

4.  Nonlinear  scattering  of  radio  waves 

Nonlinear  scattering  of  radio  waves  by  ionospheric  irregularities  was  also  investigated. 
The  details  of  the  research  were  presented  in  last  year’s  report  to  RADC.  This  process  has 
been  based  on  to  interpret  the  observed  spectral  broadening  of  VLF  radio  waves  traversing 
the  ionosphere.  In  addition,  it  has  been  shown  that  ionospheric  heating  by  HF  heater  waves 
can  be  further  enhanced  by  the  nonlinear  scattering  of  VLF  waves.  Two  papers  describing 
results  of  the  research  were  recently  published  in  technical  journals  by  Dr.  Lee  and  his 
Students.  I  hey  are  included  in  this  report  as  Appendices  B  and  C. 
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Figure  1. 


Radio  signals  transmitted  from  an  ionosonde 
have  two  different  ray  trajectories  in  the 
presence  of  ionospheric  irregularities. 
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Oblique  Heating  Experiments 

a.  Held  experiments 


D  :  Diagnostic  radio  wave 
T  :  Heater  wave 


Figure  2.  Schematic  illustration  of  oblique  ionospheric 
heating  experiments.  A  diagnostic  wave  with  a 
frequency  slightly  less  than  the  heater  wave 
frequency  is  transmitted  along  the  same  path  as 
the  heater  wave. 
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Appendix  A 


Anisotropic  Effects  of  Ionospheric  Irregularities 
on  Faraday  Polarization  Fluctuations 

Abstract 

The  previous  model  [Lee  et  al.,  1982]  of  the  Faraday  polarization  fluctuations  (FPF) 
is  extended  after  taking  into  account  the  anisotropic  nature  of  the  commonly  observed, 
rod- like  and  sheet- like  ionospheric  irregularities.  Striking  anisotropic  effects  are  seen  in  the 
longitudinal  radio  propagation.  However,  if  the  wave  propagation  angle  is  not  small  (  say 
greater  than  5®),  the  anisotropic  effects  of  rod-like  irregularities  weaken  significantly  while 
those  of  sheet-like  irregularities  remain  prominent.  In  general,  the  sheet-like  ionospheric 
irregularities  are  more  effective  than  the  rod-like  ionospheric  irregularities  in  causing  the 
FPF  of  radio  waves.  It  is  expected  that  intense  FPF  of  VHF  radio  signals  can  be  observed 
not  only  near  the  equatorial  anomaly  but  also  in  the  auroral  region. 
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1.  Introduction 


The  phenomenon  of  Faraday  Polairization  Fluctuations  (FPF)  has  been  observed  with 
low-orbit  satellite  beacon  signals  transmitted  at  the  frequencies  of  20,  40,  and  54  MHz, 
and  also  with  geostationary  satellite  signals  at  136  MHz  [Lee  et  al.,  1982  and  references 
therein).  It  results  from  the  diffractive  scattering  of  radio  waves  by  ionospheric  density 
irregularities  [Yeh  and  Liu,  1967;  Lee  et  al.,  1982;  Bhattacharyya  and  Rastogi,  1987). 
Ionospheric  irregularities  characterized  by  power-law  spectra  can  cause  the  FPF  of  radio 
signals  in  a  broad  range  of  frequencies  (  including  the  VHF  signals).  By  contrzist,  those 
with  Gaussian  spectra  can  only  give  rise  to  the  FPF  of  radio  waves  at  frequencies  up  to 
the  HF  band. 

In  the  previous  work,  the  ionospheric  irregularities  EU"e  2issumed  to  be  isotropic  for 
simplicity,  and  modelled  by  correlation  functions  having  the  same  correlation  length  in 
all  directions.  In  fact,  ionospheric  irregularities  elongate  along  the  earth’s  magnetic  held. 
They  are  highly  anisotropic  because  their  scale  lengths  adong  the  geomagnetic  field  are 
much  greater  than  those  across  the  geomagnetic  field.  Rod-like  and  sheet-like  are  the  two 
forms  of  density  irregularities  commonly  seen  in  the  ionosphere.  The  different  geometric 
effects  of  these  ionospheric  irregularities  can  be  exhibited  in  the  scintillation  measurements 
of  radio  signads  [Costa  et  al.,  1988]. 

The  purpose  of  the  present  work  is  to  investigate  the  anisotropic  effects  of  ionospheric 
irregularities  on  the  Faraday  polarization  fluctuations  (FPF)  of  radio  waves.  The  analysis 
of  FPF  presented  in  the  next  two  sections  is  done  for  both  the  rod-like  and  sheet-like 
ionospheric  irregularities.  Discussion  and  comparison  of  the  prominent  anisotropic  effects 
introduced  by  these  ionospheric  irregularities  on  the  FPF  are  subsequently  given. 

2.  Variance  of  the  Faraday  Polaurization  Fluctuations 

Consider  the  propaigation  of  a  linearly  polarized  raulio  signal  through  the  ionosphere 
modelled  ais  a  plasma  slab  of  mean  density  N  and  thickness  D.  A  Cartesiam  system  of 
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coordinates,  as  illustrated  in  Figure  1,  is  chosen  at  the  transmitter  5(0, 0,0)  with  its  z 
cixis  along  the  signal  path  from  the  transmitter  to  the  receiver  i2(0,0,  Zq)  on  the  ground. 
Ionospheric  density  irregularities,  denoted  by  A/^,  fluctuate  about  a  mean  value  <  > 

where  is  the  ratio  of  the  ambient  plasma  frequency  to  the  signal  frequency.  The  variance 
of  Faraday  polarization  fluctuations  has  been  formulated  for  single,  small-angle  forward 
scattering  cis 


(|/|^)  = 


—  /  / 
Jv  Jv 


A~  -  ^2) 

^z\Z2{Zq  -  zi){Zo  -  22) 


+  yi^)Zo 
.  2zi{Zo  -  21 ) 


k{x2^  +  y2^)Zo'\  _.^reA!(xi^  +  y\^)ZQ-\  _.^r£A:(x2^  +  y2^)^0~ 

2z2{Zq-Z2)  J  L  2zi{Zq-zi)  I  2z2{Zq  -  Z2)  ■ 


(1) 


where  C(ri  —  T2)  =  ((^/^(n))(^7^(^2))*)  is  the  correlation  function  of  the  density 
irregularities  (see  equation  (19)  in  Lee  et  al.  [1982]). 


3.  Rod-like  ionospheric  irregularities 

The  anisotropy  of  ionospheric  irregularities  cam  be  accounted  for  by  a  longer  correlation 
length  in  the  direction  parallel  to  the  m^netic  field  and  a  shorter  correlation  length  in 
the  perpendicular  direction.  Ionospheric  irregularities  with  a  rod-like  configuration  can  be 
reaisonably  expressed  by  the  following  correlation  function 


C(ri  -  r2)  =  {Wn?) 


where  ((A/^)^)  corresponds  to  the  mean  squared  density  fluctuations,  ^j.(^||)  is  the 
correlation  length  across  (along)  the  eamth’s  magnetic  field.  The  system  of  the  prime 
coordinates  is  designated  in  Figure  2. 

Since  equation  (1)  is  written  in  the  unprime  coordinates,  the  correlation  function  (2) 
needs  also  be  expressed  in  the  same  coordinates.  The  prime  and  unprime  coordinates  are 
related  by 

sin  V'  —  cos  V'  0  \  (  ^\ 

cos  d  cos  cos  9  sin  V*  —  sind  1  I  1/  1 

sin  6  cos  V'  sin  9  sin  ^  cos  9  /  \  z  / 
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It  is  convenient  to  make  the  following  coordinate  transformation  to  carry  out  the 


integrations  in  equation  (1)  ; 


<  =  II  -  12 


^~y\-y2 


C  =  21  -  22 


2X  =  ii  -h  12  2Y  ~  y\  +y2  2Z  z\  ->r  zo 
After  the  above  transformation,  the  correlation  function  has  the  form  of 

C{i,-nX)  =((A7^)^)exp  ^(4sinV>  -  qcos^j 


exp  —  |4cosdcosi/>  +  qcosdsini/)  —  (^sin^l 


•  exp  —  sin  &  cos  ij)  +t)  sin  d  sin  4>  +  cos  0\ 

[h 

Integrations  of  equation  (1)  over  x  and  y  with  infinite  limits  yield 

Jb2  r-a+Dsecx 


/Dxcx  /*oo  /*oo 

d(  dr,  diCii,-n,0  (6) 

■  Dsccx  J—oo  J —CO 

■  X;G(«)  sin  [G(i)(«2  +  -,2)] 
t 

where  G{t)  zind  are  defined  as  follows 

G{t)  =  Zok  [2C(2Z  -  Zo)  -h4f(ZZo  -  Z2  -  C^/4)]"^  (7) 

Y,  h{t)  =  2h{0)  -  h{e)  -  h{-t)  (8) 

t 

To  simplify  the  integrations  in  equation  (6),  ($,  q,  C)  is  transformed  into  (u,  u,  w).  They 
are  related  by  the  same  trzmsformation  tensor  in  equation  (3).  This  transformation  converts 


equation  (6)  into 


(Kr)  = 


16.(1-  <  f%  >)^  Ja 


ra+Dsecx 


dZ  /  dv 


■  /I  ("?)  /I 

■  ^  |<?(f)|  sin  pG(t)|((u  COS0  -b  ti;sin0)^  -I-  u^)j 
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where  Va  =  wcotd-  D  sec  x  esc  6,  vf,  =wcot6  +  D  sec  x  esc  9,  and  all  the  C’s  in  (7)  have 
been  replaced  by  (— u  sin 9  +  w  cos  9). 


To  cairry  out  the  integrations  in  equation  (9)  further,  the  following  treinsformations  are 
made: 


V  =  V  cos  9 


W  =  w  sin  9 


(10) 

V  -\-W  =  3  V  -W  =  p 

This  set  of  trcinsformations  renders  the  integration  over  p  trivieil.  The  result  after  the 


integration  over  p  is 


O  k'^UAflf)  r 
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sin0 


exp  - 


cos  9 


£ii  sin  9 


exp 


^|G(()|sin  [|G(1)|(«2  +  ,2)] 


The  integrations  over  u  and  a  can  be  carried  out  exactly.  Then,  the  integration  over  Z  can 
simply  be  approximated  by  the  product  of  the  interval  of  integration  and  the  integrand 
wherein  all  the  Z’s  axe  set  equal  to  Zm,  the  distance  between  the  tramsmitter  and  the 
mean  height  of  the  plaisma  slab. 

Finally,  the  variance  of  the  Faraday  rotation  fluctuations  is  obtained  as 

<  -fl  >2 


<|/r)  = 


4(1-  <  >f 


{{AN/NY) 


Dsecx 


sin^  9  cos^  9 


-1 


(12) 


where 


sind 


[1  -  2Fs{L)Fc{Ls)  -  2Fc{L)F,{L,)] 


cos  9 
in 


[1  -  2Fs{L)FciLc)  -  2Fc{L)Fs{Lc)] 


(13) 


in  which  ,  lY  =  = 


cd 


,  d  is  the  mean  height  of  the 


density  irregularities,  and  the  functions  Fj(®)  and  Fc(x)  are  defined  by  equations  (14)  and 
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(15),  respectively. 


Fa{x)  =  i[cos(x^)  +  sin(x~)]  —  [cos(i^)C(x)  +  sin(x^)S{i)]  (14) 

Fc{x)  =  ^[cos(x^)  —  sin(x^)]  —  [cos(x^)5(x)  —  sin(x^)C(x)]  (15) 

where  5(x)  and  C{x)  are  the  Fresnel  sine  and  the  Fresnel  cosine  functions,  respectively. 

The  expression  for  FPF  variajice  in  (12)  can  be  very  much  simplified  because  the 
arguments  of  the  functions  F3{x)  and  Fc{x)  are  very  small  with  the  adopted  ionospheric 
parameters  discussed  in  the  next  section.  Therefore,  Fj(x)  and  Fc(x)  can  be  expanded 


into  power  series,  and  terms  of  the  order  of  x  a^e  kept  to  yield  the  result 

0  ^  <  /w  (ekdV^^D  secx  o 

{|/|2)  = - .7-  '  ■: - ^ - {{AN/Nf)  (16) 

5(1-  <  >)2[sin6  +  (^i/^llj  cos6] 

which  is  valid  for  in  sind  ^  yj{ed/k)  amd  cosB  \/{tdfk).  The  term  associated  with 


the  (fi/^|j)  in  (16)  accounts  for  the  anisotropic  effects  of  rod-like  irregulauities  on  the  FPF. 
In  the  case  of  exactly  longitudinal  propagation,  the  FPF  variance  cam  be  calculated  from 
(11)  by  setting  d  =  0®.  For  £x  cosO  »  i/(cd/fc),  the  result  is 


2  _  k<  fjf  >^{ekd)^^^D  sec  X 


(1/r)  = 


5(1-  <  >)2 


(«||/«j.)  ((A.V/JV)2) 


4.  Sheet-like  ionospheric  irregularities 

Sheet-like  ionospheric  irregularities  au-e  seen  in  the  auroral  oval,  while  rod-like 
irregularities  at  higher  and  lower  latitudes  [e.g.,  Aau’ons,  1982].  Sheet-like  irregularities 
cam  be  represented  by  the  correlation  function  in  the  form  of 


These  sheets  may  coincide  with  the  meridian  plane  or  a  plane  which  is  orthogonad  to  the 
meridiam  plame.  In  the  former  case,  the  y' z*  plame  represents  the  meridian  plame  amd  the 
x'  aixis  points  to  the  east-  west  direction.  Whereas,  in  the  latter  case,  the  x'  axis  is  along 
the  north-south  direction  and  the  y'z'  plane  intersects  the  meridian  plane  at  right  angles. 
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The  corresponding  variamces  of  the  FPF  caused  by  the  sheet-like  irregulcirities  are 


found  to  be 


2  fc  <  /^  >^  (eA;d)^/^Dsecx  /AiVX^ 


<  i/r  >= 


5(1—  <  /^  >)2(sin0  -f  cos^) 
and,  for  the  longitudinal  propagation  {$  =  0°), 


(?)<“)■ 


2  ^  <  fit  (tA:d)^/^Dsecx 


<  l/f  >~ 


5(1-<7J>)2 


(«||/£j.)  <  {AN /NY 


where  £||  »  is  assumed.  Compare  equations  (19)  aind  (20)  with  (16)  and  (17),  one  can 
see  the  distinctive  anisotropic  effects  introducd  by  the  sheet-like  and  rod-like  ionospheric 
irregularities  on  the  FPF  of  radio  waves.  It  is  evident  from  (16)  that,  as  sin  9  >  (£j./^||)  cos  0 
[e.g.,  9  >  5.7°  for  ^x/^||  =  0.1],  the  anisotropic  effect  of  rod-like  ionospheric  irregularities 
is  no  longer  prominent  for  large  propagation  angles.  By  contrast,  the  anisotropic  effect  of 
sheet-like  irregularities  does  not  depend  upon  the  propagation  angles  strongly  zu:cording 
to  (19).  The  most  prominent  anisotropic  effects  of  ionospheric  irregularities  are  seen  in 
the  longitudinal  radio  propagation.  Further,  both  the  rod-like  and  sheet-like  irregularities 
exhibit  the  same  effects  on  the  longitudinal  wave  propagation. 


5.  Discussion  and  Conclusion 

We  have  extended  the  previous  model  of  the  Faraday  Poleuization  Fluctuations  (FPF) 
by  taking  into  account  the  anisotropy  of  ionospheric  irregularities.  The  anisotropic  effects 
of  the  commonly  observed,  rod- like  and  sheet-like  ionospheric  irregulau'ities  on  the  FPF 
of  radio  waves  are  exzunined.  The  striking  anisotropic  effects  of  both  types  of  ionospheric 
irregularities  are  expected  for  the  longitudinal  radio  propagation. 

However,  if  the  wave  propagation  angle  (d)  is  not  small  (say,  9  >  5°),  the  anisotropic 
effects  of  rod-like  irregularities  weaken  significantly,  while  those  of  sheet-like  irregularities 
remain  prominent.  When  (£^/£||)co8d  is  much  less  than  sind,  equation  (16)  reduces  to  a 
form  similar  to  equation  (27)  of  Lee  et  al.  [1982].  Hence,  the  refractive  wave  scattering  by 
rod-like  ionospheric  irregularities  for  a  large  prop£^;ation  angle  can  be  adequately  described 
by  the  isotropic  random-  medium  model  of  Lee  et  al.  [1982]. 
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The  primary  distinction  between  equation  (27)  of  Lee  et  al.  [1982]  and  equation 
(19)  is  the  large  factor,  (^|j/fj.)  ~  0(1)  which  cheiracterizes  the  anisotropic  effects  of 
sheet-like  irregularities.  Because  of  this  large  factor,  the  ionospheric  conditions  for  the 
occurrence  of  FPF  cam  be  signiAcantly  relaxed.  For  instance,  the  FPF  caused  by  the 
sheet-like  ionospheric  irregularities  requires  less  density  fluctuation,  and/or  background 
plasma  density  to  occur.  We,  therefore,  conclude  that  the  sheet-like  irregularities  axe  more 
effective  tham  the  rod-like  irregularities  in  causing  the  FPF  of  rauiio  waves.  Since  the  sheet¬ 
like  irregulau'ities  appeair  in  the  auroral  oval  while  the  rod-like  irregularities  at  higher  and 
lower  latitudes  [Aarons,  1982],  the  present  work  predicts  that  intense  FPF  of  VHF  signals 
cam  be  observed  not  only  near  the  equatoriad  anomady  but  adso  in  the  auroral  region. 
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Figure  1. 


Figure  2. 


Figure  Captions 


Geometry  of  the  tramsionospheric  radio  propaigation  in  the  presence  of 
ionospheric  density  irregularities. 


Illustration  of  ainisotropic  ionospheric  density  irregularities. 
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Appendix  B 

Specticd  Broadening  of  VLF  Radio  Signals 
Traversing  the  Ionosphere 
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Nearly  monochromatic  lignalt  at  13.6  kHx  ±1  Hz  irvjected  from  a  ground-baaed  VLP  tranamitter 
can  experience  a  bandwidth  expanaion  aa  high  aa  1%  (-s-  100  Hz)  of  the  incident  wave  frequency 
aa  they  traveric  the  ionoaphere  and  reach  aateliite  altitudea  in  the  range  of  600-3800  km  (Bell 
et  al.,  1983)  The  off-carrier  componenta,  having  electroatatic  nature,  are  believed  to  be  induced 
lower  hybrid  wave  modea.  We  inveatigate  two  different  aource  mechaniama  that  can  potentially 
result  in  the  obaerved  spectral  broadening  of  ii\iected  monochromatic  VLF  waves.  One  is  the 
nonlinear  scattering  of  VLF  signals  by  ionospheric  density  fluctuations  that  renders  the  nonlinear 
mode  conversion  of  VLF  waves  into  lower  hybrid  waves.  These  electrostatic  modes  result  when 
the  iiviected  VLF  waves  are  scattered  by  ionospheric  density  fluctuations  with  scale  lengths  less 
than  V^w(c/ui,<)(D./u>o)^^^  2  km  in  the  upper  ionosphere,  where  c,  uiy«,  D.,  and  Us  are  the 

speed  of  light  in  vacuum,  the  plasma  frequency,  the  electron  cyclotron  frequency,  and  the  VLF 
wave  frequency,  respectively.  In  the  absence  of  ionospheric  irregularities,  a  second  mechanism  that 
involves  a  parametric  iiutabilicy  can  excite  the  lower  hybrid  waves  (Lee  and  Kuo,  1984).  This 
process  tends  to  produce  a  spectrally  broadened  transmitted  pulse  with  peaks  at  a  discrete  set  of 
frequencies  on  both  tides  of  the  nominal  carrier  frequency.  By  contrast,  the  nonlinear  scatteriirg 
mechanism  generates  single-peaked  spectra  centered  at  the  carrier  frequency.  Both  types  of  spec¬ 
tra  were  obaerved  in  the  experiments.  Therefore,  the  two  suggested  source  mechanisms  coirtnbute 
additively  to  the  observed  spectral  broadening  of  it\iected  VLF  waves. 


1.  Introduction 

Nearly  monochromatic  signals  at  13.6  kHz  ±  1  Hz  in¬ 
jected  from  a  ground-based  VLF  transmitter  were  recently 
observed  to  experience  bandwidth  expansion  as  they  tra¬ 
versed  the  ionosphere  and  reached  the  satellite  altitudes 
in  the  range  of  600-3800  km  [Bell  et  al.,  1983].  This  ex¬ 
pansion  of  bandwidth,  which  results  in  a  proportional  in¬ 
crease  in  signal-to-noise  ratio,  may  be  aa  large  aa  1%  (— 
100  Hz)  of  the  carrier  frequency,  and  the  off-carrier  compo¬ 
nents  are  thought  to  be  electrostatic  in  nature  [/non  and 
Bell,  1985].  They  occur  only  in  the  presence  of  impul¬ 
sive  VLF  hiss  and/or  a  lower  hybrid  resonance  (LHR)  noise 
band  with  an  irregulau:  cutoff  frequency,  and  only  for  signals 
whose  frequencies  exceed  the  LHR  frequency  at  the  satel¬ 
lite  location.  Based  on  these  observed  characteristics  of  the 
spectrally  broadened  transmitter  signals,  one  may  suspect 
the  off-carrier  components  to  be  either  electroatatic  whistler 
wave  modea  or  lower  hybrid  wave  modes. 

The  linear  scattering  source  mechanism  producing  quaai- 
electrostatic  whistler  wave  modes  was  first  suggested  by 
Belt  et  al.  [1983],  who  hypothesized  the  creation  of  the 
required  ionospheric  density  fluctuations  by  precipitating, 
low-energy  ( <  1  keV)  electrons.  Such  precipitation  events 
have  been  accompanied  by  both  VLF  hiss  and  irregular 
LHR  noise  bands  [McEuien  and  Barrington,  1967;  Laaspere 
et  al.,  1971;  Gumett  and  Prank,  1972].  Bell  et  oL  [1983] 
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then  speculate  that  the  broadening  of  the  transmitted  pulse 
spectrum  results  from  the  scattering  of  the  initial  signals 
from  the  precipitation-induced  density  fluctuations  and  the 
subsequent  coupling  into  quasi-electrostatic  whistler  wave 
modes  of  short  wavelength.  The  Doppler  shift  associated 
with  these  short-wavelength  modes  is  estimated  to  be  large 
enough  to  produce  the  bandwidth  expansion  of  the  signals 
measured  on  a  moving  satellite. 

However,  we  interpret  the  off-carrier  components  of  the 
spectrally  broadened  transmitter  signals  to  be  lower  hybrid 
wave  modes.  The  main  reason  is  that  the  lower  hybrid 
wave  modes  (eigenmodes  of  ionospheric  plasmas)  experience 
less  spatial  attenuation  than  the  quasi-electrostatic  whistler 
wave  modes  (quasi-modes  of  ionospheric  plasmas)  during 
propagation  from  the  source  regions  to  the  satellite  loca¬ 
tion.  Two  source  mechanisms  are  proposed,  whereby  the 
transmitted  VLF  waves  can  potentially  produce  the  lower 
hybrid  waves.  One  mechanism  comparable  to  that  of  Bell  et 
al.  [1983],  termed  “nonlinear  scattering  process,”  occurs  in 
the  presence  of  ionoepheric  density  irregularities.  The  other 
one,  by  contrast,  does  not  require  the  presence  of  ionoepheric 
irregularities.  It  is  a  parametric  instability  that,  excited  by 
intense  enough  VLF  waves,  can  generate  lower  hybrid  waves 
and  ionospheric  irregularities  concomitantly  [Lee  and  Kuo, 
1984]. 

The  first  mechanism,  presented  in  section  2,  differs  from 
that  of  Bell  et  al.  The  newly  proposed  nonlinear  process  is  a 
scattering  of  the  whistler  mode  VLF  signals  by  preexisting 
density  fluctuations  that  render  a  mode  conversion  to  lower 
hybrid  waves,  instead  of  a  linear  process  that  couples 
energy  into  the  quasi-electrostatic  branch  of  the  whistler 
mode.  We  formulate  a  theory  of  VLF  wave  scattering  off 
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Ionospheric  irregularities,  showing  how  this  process  depends 
upon  the  local  fractional  density  fluctuations,  the  ratio 
of  their  scale  lengths  to  the  VLF  wavelength,  etc.  The 
scattering  of  VTF  waves  by  ionospheric  density  fluctuations 
generally  causes  elllptically  polarized  modes.  The  induced 
elliptically  polarized  modes,  however,  may  be  predominantly 
electrostatic  under  certain  ionospheric  conditions. 

The  second  mechanism,  described  in  section  3,  produces 
a  spectrally  broadened  signal  via  a  parametric  instability 
[Lee  and  Kuo,  1984]  The  injected  VLF  wave,  if  intense 
enough,  parametrically  excites  a  magnetic  field-aligned 
purely  growing  mode  together  with  two  sidebands  of  lower 
hybrid  waves.  The  parametric  excitation  process  requires 
the  frequency  and  wave  vector  matching  relations  to  be 
satisfied.  The  two  lower  hybrid  wave  sidebands  are 
characterized  by  their  wave  vectors  which  are  the  sum  and 
difference  of  the  wave  vectors  of  the  whistler  pump  and  the 
purely  growing  mode,  respectively. 

The  lower  hybrid  waves  generated  by  either  of  the  two 
proposed  mechanisms  can  have  wavelengths  much  shorter 
than  that  of  the  incident  VLF  waves.  In  other  words,  much 
greater  wave  numbers  are  introduced  by  the  suggested  source 
mechanisms  for  the  injected  VLF  waves  that,  essentially, 
become  electrostatic  modes.  These  electrostatic  plasma 
modes  with  large  wave  numbers,  detected  by  satellites 
moving  across  the  geomagnetic  field,  give  rise  to  Doppler 
shifts  that  were  observed  as  the  broadened  spectra  (i.e., 
off-carrier  components)  of  the  injected  VLF  waves.  The 
preliminary  results  were  presented  in  Lee  et  al.  [1986]. 

2.  Nonlinear  VV ave  Scattering 


Theory 

A  monochromatic  V'LF  wave  transmitted  from  a  ground- 
based  station  into  space  has  been  observed  to  change  from 
linear  into  circular  polarization  (i.e.,  whistler  mode).  If 
a  ducted  whistler  wave  mode,  propagating  exactly  along 
the  geomagnetic  field,  is  considered  for  simplicity,  the  wave 
electric  field  may  be  represented  as 

E  =  Eolx -h  iyjexp[:(leo2  -  Wot)]  (1) 

where  the  z  axis  has  been  taken  along  the  geomagnetic  field 
lines,  IS  the  transmitted  wave  frequency,  and  ko  the 
associated  wave  vector,  assumed  to  be  along  the  z  axis  as 
illustrated  in  Figure  la.  Propagating  into  the  ionosphere 
with  no  irregular  structure,  the  whistler  wave  satisfies  the 
electromagnetic  wave  equation 


O 


c2 


(2) 


where  jo  is  the  uniform  oscillatory  current  driven  by  the 
incident  wave  field  in  the  ionospheric  plasma. 


le^NoEo 
m(LV„  -  n.  ) 


and  e.  No,  m,  and  Q,  are  the  electron  charge,  uniform 
background  plasma  density,  electron  mass,  and  the  unsigned 
electron  cyclotron  frequency,  respectively.  In  the  presence  of 
field-aligned  ionospheric  density  irregularities,  the  scattered 
wave,  allowing  for  electrostatic  modes,  can  be  described  by 


V(V  E.)  -  V^E. 


c2  dt^^‘ 


(3, 


where  ij,  is  the  nonuniform  current  induced  by  the 
interaction  of  the  whistler  wave  field  with  the  density 
irregularities.  We  take  the  i  axis  along  the  direction  of  the 
density  irregularities  which  are  assumed  to  be  of  the  form 


Sn  =  Sn  exp[ifcij 


( 4 ) 


depicted  graphically  in  Figure  16.  The  nonuniform  current 
6j,  can  then  be  expressed  as 


i5j.  =  -e(NoSv  ■+■  ^jo) 

fy  o 


( 5  / 


where  6v  is  the  induced  velocity  perturbation.  The  scattered 
wave  field  has  the  general  form  of 


E,  =  [xE,  -e  tyEyjexplt(koZ  -  uJot)]  1 6) 


Solving  the  wave  equation  together  with  the  electron 
momentum  equation,  the  scattered  wave  field  is  found  to 
be  elliptically  polarized, 

E.  =  Eo(Sn/No)[x  +  (x  -t  iy)Q]exp[t(ftoZ  -  ui„£)]  (7) 


where  a  is  proportional  to  the  amplitude  of  the  circularly 
polarized  component  of  the  wave  defined  by 


Q 


n.) 


(8) 


where  lOp,  and  c  are  the  electron  plasma  frequency  and 
the  speed  of  light  in  vacuum,  respectively,  and  k  is  the 
wave  number  of  the  density  perturbation  defined  by  (4) 
Rewriting  the  scattered  field  in  terms  of  the  ionospheric 
irregularity  scale  length,  \  =  2rr/fc,  we  find 


Fig.  1  Geometry  of  wave  propagation  for  (a)  uniform 
lonoaphere  and  (b)  ionosphere  with  irregularities 
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E,  =  Eo  exp[ifcij(6n/fV„)  exp[n'*:o2  —  uiotlj 

>,2 

X  [^_x_^-r(x  +  zy)— ■  (9) 

LP  CP 

where 

Ac  =  y27r(c/ui,.)(D./u;„)‘'^  5:  4.4.'fc„  (10) 

The  scattered  field  is  composed  of  a  linearly  polarized  (LP) 
part  and  a  circularly  polarized  (CP)  part  as  indicated  in  (9). 
If  the  ratio  (A^/Aj)  <  1  (i.e.,  |k|  >  |ko|),  the  scattered  wave 
is  dominated  by  the  linearly  polarized  component  which, 
corresponding  to  a  lower  hybrid  mode,  oscillates  in  the 
direction  of  the  irregularity  wave  vector  k;  hence,  the  wave 
is  electrostatic  in  nature  and  exhibits  a  broadened,  enlarged 
wave  vector  spectrum  relative  to  the  incident  wave.  The 
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whistlers  are  not  eigenmodes  but  quasi-modes  of  ionospheric 
plasmas.  Therefore,  during  the  ionospheric  propagation 
of  quasi-electrostatic  whistlers,  the  ionospheric  plasmas  do 
not  respond  in  phase  to  these  quasi-modes.  This  phase 
mixing  (or  phase  mismatching]  can  attenuate  the  waves 
significantly  during  their  propagation  from  the  source  region 
to  the  satellite  location. 

By  contrast,  our  proposed  nonlinear  scattering  mechanism 
can  produce  lower  hybrid  waves  which  are  eigenmodes  of 
ionospheric  plasmas,  experiencing  less  spatial  attenuation 
than  the  quasi-electrostatic  whistler  modes  during  propaga¬ 
tion.  This  mechanism  is  effective  provided  that  field-aligned 
ionospheric  irregularities  exist  and  that  the  VLF  wave  fre¬ 
quency  is  greater  than  the  local  LHR  (lower  hybrid  reso¬ 
nance)  frequency.  The  characteristics  of  the  nonlinear  scat¬ 
tering  mechanism,  then,  explain  the  observed  facts  that  the 
spectral  broadening  phenomenon  occurs  only  in  the  pres¬ 
ence  of  impulsive  VLF  hiss  and/or  a  LHR  noise  band  with 
Fig.  2.  Sdiematic  of  Doppler  shift  phenomena.  an  irregular  cutoff  frequency,  and  only  for  signals  whose  fre¬ 

quencies  exceed  the  LHR  frequency  at  the  satellite  location. 


Doppler  shift  frequency  measured  by  a  moving  satellite  will 
be  enlarged  as  elaborated  in  the  next  section. 

Anaiyji) 

The  condition  found  above  for  effective  nonlinear  scatter¬ 
ing  into  electrostatic  modes  requires  A  A,.  Typical  plasma 
parameters  for  the  upper  ionosphere/lower  magnetosphere 
region  of  interest  are  w,,/2ir  =  0.7  MHz  and  fl,/2ir  =  0.6 
MHz.  Assuming  an  incident  frequency  of  a/o/2jr  =  13.6  kHz, 
we  calculate  Ae  %  2  km.  The  magnitude  of  the  bandwidth 
expansion  of  the  scattered  VLF  wave  measured  by  a  moving 
satellite  due  to  the  Doppler  shift  effect  is  given  by 

k.  v,  v.coaff  V, 

where  v,  is  the  velocity  of  a  satellite  moving  across  the 
Earth’s  magnetic  field  as  shown  in  Figure  2,  and  k,  = 
kx  kai.  Given  v,  ~  8  km/s  and  A/  ~  100  Hz,  we  find  A  ~ 
80  m  Ac  ~  2  km,  so  the  condition  for  nonlinear  scattering 
is  well  satisfled.  Thus  we  have  found  that  ionospheric 
irregularities  with  scale  lengths  of  several  tens  of  meters 
can  produce  Doppler  shifts  of  the  observed  magnitude. 
Since  fc(A)  3>  (<^)ko(^o),  the  simplification  made  in  the 
formulation  of  the  theory,  namely,  the  assumption  that 
the  ducted  whistlers  propagate  exactly  along  the  Earth’s 
magnetic  field,  is  justified.  Under  this  condition,  even  if 
the  ducted  whistlers  have  nonzero  wave  normal  angles,  the 
formulation  of  nonlinear  scattering  mechanism  will  not  be 
significantly  modified. 

Ionospheric  irregularities  can  occur  naturally.  There 
are  many  mechanisms  that  can  generate  field-aligned  iono 
spheric  irregularities  [e  g.,  fejer  and  Kelley,  1980],  Particle 
precipitation  invoked  by  Bell  et  al.  [1983]  can  produce 
irregularities  with  a  broad  range  of  sesde  lengths  in  a  region 
near  the  F  layer  peak,  about  300  km  above  the  Earth’s 
surface.  The  F  region  irregularities,  as  Bell  et  al.  suggested, 
provide  a  favorable  ionospheric  condition  for  the  linear 
scattering  of  the  VLF  waves  and  the  subsequent  coupling 
into  quasi-electrostatic  whistler  wave  modes.  However,  the 
amplitude  of  such  a  quasi-mode  wave  decreases  with  distance 
as  it  propagates  upward  to  the  satellite  altitudes  of  600- 
3800  km.  This  arises  from  the  fact  that  the  electrostatic 


3.  Parametric  Excitation  op  Lower 
Hybrid  Waves 

In  this  section  we  discuss  the  proposed  second  mechanism 
responsible  for  the  observed  VLF  spectral  broadening. 
While  the  nonlinear  scattering  mechanism  does  not  work  in 
the  absence  of  ionospheric  irregularities,  lower  hybrid  wave 
modes  can  be  excited  parametrically  by  the  VLF  waves  in 
the  ionosphere.  The  parametric  excitation  of  lower  hybrid 
waves  by  the  injected  VLF  waves  has  been  suggested  by  us 
[lee  and  Kuo,  1984]  with  the  following  scenario. 

Propagating  from  the  neutral  atmosphere  into  the  iono¬ 
sphere,  the  incident  VLF  wave  can  be  considered  to  be  a 
ducted  whistler  mode  with  ko  =  k^z  along  the  geomagnetic 
field  lines,  taken  to  be  the  z  axis  of  a  Cartesian  system  of  co¬ 
ordinates.  Assuming  the  space-time  dependence  of  the  per¬ 
turbations  to  be  exp[i(k  •  r  -  wt)],  the  proposed  parametric 
instability  can  be  described  by  the  following  wave  frequency 
and  wave  vector  matching  relations: 

Wx.f  —  Wf  —  Wo  ~  wi  —  "4-  Wj  (12o) 

ki4  -  xk,  =  zko  =  ki-  +  xk,  {12b) 

where  the  subscripts  minus/plus  and  s  represent  the 
Stokes/anti-Stokes  components  of  the  high-frequency  lower 
hybrid  sidebands  and  the  field-aligned  zero-frequency  mode, 
respectively.  By  the  definition  of  (12a),  the  Stokes  and 
anti-Stokes  components  refer  to  the  lower  hybrid  sidebands 
whose  wave  frequencies  are  downshifted  and  upshifted, 
respectively,  by  w,  with  respect  to  the  whistler  wave 
frequency  (wo)-  For  a  purely  growing  mode  (i.e.,  Re(w,)=0), 
the  Stokes  and  anti-Stokes  components  of  lower  hybrid  waves 
are  characterized  by  their  wave  vectors,  shown  in  (12b), 
which  are  the  difference  and  sum  of  the  wave  vectors  of 
the  whistler  and  the  purely  growing  mode,  respectively. 
Both  the  Stokes  and  the  anti-Stokes  components  of  lower 
hybrid  waves  have  to  be  considered  in  this  two-dimensional 
instability  process,  since  |k,|  ^  |ko|  and  the  concomitantly 
excited  low-frequency  mode  is  a  purely  growing  mode.  Thus, 
jki+l  ~  |ki_|  ~  |k.|;  it  is  found  that  the  lower  hybrid  waves 
have  wave  vectors  large  enough  to  yield  the  observed  spectral 
broadening  due  to  the  Doppler  effect. 

It  should  be  pointed  out  that  a  four-wave  rather  than  a 
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simpler  three-wave  interaction  process  [Berger  and  Perkins, 
1976;  Riggin  and  Kelley,  1982]  is  proposed  by  Lee  and 
Kuo  [1984].  While  the  proposed  parametric  instability  has 
a  broad  range  of  VLF  (whistler)  wave  frequency  Wo,  the 
three- wave  interaction  process  [Berger  and  Perkins,  1976]  is 
restricted  to  a  very  narrow  frequency  band,  namely,  < 
uJo  <  where  u)^ff  is  the  lower  hybrid  resonance 

frequency  Further,  the  excited  low-frequency  mode  in  the 
work  of  Lee  and  Kuo  [1984]  is  a  zero-frequency  mode.  This 
four-wave  interaction  mechanism  is  found  to  require  much 
lower  thresholds  than  the  three-wave  interaction  mechanism 
to  excite  lower  hybrid  waves  in  the  ionosphere.  Based  on  this 
proposed  four-wave  interaction  process,  we  can  interpret  the 
observations  of  the  VLF  wave-excited  lower  hybrid  waves  in 
the  ionosphere  indicated  in  the  experiments  of  the  Franco- 
Soviet  ARCAD  3  satellite  [Berthelier  et  al.,  1982]  and  others 
[e  g.,  Chmyrev  et  al.,  1976]. 

.At  this  point,  we  should  comment  on  the  effectiveness 
of  the  parametric  instability  in  causing  the  observed  VLF 
spectral  broadening  phenomenon.  Although  the  production 
r.f  lower  hybrid  waves  by  the  parametric  instability  does 
not  require  the  presence  of  naturally  occurring  ionospheric 
irregularities,  this  mechanism  needs  a  threshold  to  operate. 
The  thresholds  in  terms  of  whistler  (VLF)  wave  electric  fields 
can  be  calculated  from  equation  (22b)  of  Lee  and  Kuo  [1984]: 


=  0.86(fcJuf,/fl,)(m/e)|rj| 


1  4-  1  -h 


1/2 


(13) 


where  k,  v,,,  u,,  and  fl.  are  the  wave  number  of  field- 
aligned  ionospheric  irregularities,  the  electron  thermal 
velocity,  the  electron-ion  collision  frequency,  and  the  electron 
cyclotron  frequency;  rj  is  defined  by  [\  +  (M/m){ko/k)^]/[\  — 
[M lm){ko/k)^{Q,/u)f,)^]  where  M/m,  ko,  and  Wp,  are 
the  ratio  of  oxygen  ion  (O’*")  mass  to  electron  mass, 
the  wave  number  of  the  incident  whistler  wave,  and  the 
electron  pleisma  frequency,  respectively.  It  should  be  noted 
that  the  parametric  wave  number  k  is  determined  by 
the  dispersion  relation  for  the  lower  hybrid  waves,  Wo  = 
-ii,„[l  T  (M/m)(ko/k,)^Y^^,  where  is  the  lower  hybrid 
resonance  frequency  defined  by  =  wp,/(i  -f- 
If  we  choose  the  following  model  parameters;  v,,  =  1.3  x  10^ 
(m/s),  u/p,/27r  =  0.7  MHz  (i.e..  No  =  6  x  10®  m'^),  i/,  =  10 
Hz,  Q,/27r  =  0.6  MHz,  the  threshold  field  £„  is  found  to  be 
18  pV/m.  According  to  Inan  and  Bell  [1985],  the  measured 
power  flux  of  transmitted  VLF  signals  is  2  x  10“'^  W /m®  at 
the  satellite  location;  namely,  the  measured  field  intensity 
of  whistlers  is  12  pV/m,  which  is  of  the  same  order  of 
magnitude  as  the  threshold  of  the  parametric  instability. 
The  measured  power  flux  (field  intensity)  of  whistler  (VLF) 
waves  in  the  lower  ionosphere  is  2  x  10”®  W/m^  (1.2 
mV  /m),  which  is  also  of  the  same  order  of  magnitude  as 
the  instability  thresholds  calculated  by  Lee  and  Kuo  [1984]. 
Therefore,  it  is  possible  for  the  instability  to  be  excited 
though  the  VLF  transmitters  are  operated  in  pulsed  mode 
with  a  typical  duration  of  seconds. 

It  is  interesting  to  note  that  the  proposed  two  source 
mechanisms  are  expected  to  produce  broadened  VLF  spectra 
with  different  shapes.  The  parametric  instability  will 
preferentially  excite  lower  hybrid  waves  with  the  scale 
length  (A  =  27r/fc)  of  ~  90  m  at  the  satellite  location. 
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This  likely  generates  a  two-humped  spectrum  suppressed 
at  the  carrier  frequency  (/»)  relative  to  the  sidebands 
peaking  at  approximately  /o±100  Hz,  whereas  the  nonlinear 
scattering  mechanism  simply  broadens  the  bandwidth  of  the 
injected  VLF  waves,  producing  single-peaked  spectra.  This 
can  be  clearly  seen  from  equation  (9),  showing  that  the 
scattered  field  intensity  E,  is  the  incident  field  intensity 
Eo  multiplied  by  the  fractional  plasma  density  fluctuation 
(Sn/No  <  1).  This  naturally  leads  to  the  broadened 

V'LF  spectrum  peaking  at  the  carrier  frequency.  Tnese 
two  types  of  broadened  VLF  spectra  have  indeed  been 
observed.  Nevertheless,  the  single-peaked  spectra  were 
much  more  frequently  seen  than  the  double-humped  spectra 
(H.  G.  James,  personal  communications,  1988).  The 
characteristics  of  our  proposed  two  source  mechanisms  agree 
well  with  the  observations.  Hence,  the  observations  of  VLF 
spectral  broadening  offer  additional  evidence  supporting  the 
excitation  of  lower  hybrid  waves  by  V'LF  waves  in  the 
ionosphere  as  suggested  by  Lee  and  Kuo  [1984]. 

4.  Summary  and  Conclusion 

In  summary,  we  have  investigated  two  possible  source 
mechanisms  that  may  be  responsible  for  the  observed 
spectral  broadening  of  injected  V'LF  waves.  In  the  presence 
of  ionospheric  irregularities  with  various  scale  lengths  as 
short  as  several  tens  of  meters,  the  nonlinear  scattering  of 
the  VLF  waves  off  these  density  irregularities  can  produce 
electrostatic  modes  with  larger  wave  vectors  which  give  rise 
to  the  apparent  spectral  broadening  through  the  Doppler 
shift  observed  by  the  moving  satellites.  The  amplitude  of 
the  electrostatic  modes  depends  linearly  on  the  amplitude 
of  the  ionospheric  irregularities.  The  broadening  produced 
through  this  mechanism  is  of  the  same  magnitude  as 
the  observed  values.  This  mechanism  produces  single- 
peaked  broadened  spectra  of  VLF  waves,  which  were  most 
frequently  observed  in  the  experiments.  There  are  many 
natural  sources  that  can  generate  ionospheric  irregularities 
with  a  broad  range  of  scale  lengths,  providing  favorable 
ionospheric  conditions  for  the  operation  of  the  nonlinear 
scattering  process. 

In  the  absence  of  ionospheric  irregularities,  the  spectral 
broadening  of  the  incident  wave  packet  may  also  stem 
from  the  second  mechanism,  which  involves  the  parametric 
excitation  of  lower  hybrid  waves  [Lee  and  Kuo,  1984].  This 
mechanism  requires  a  threshold  that  can  be  exceeded  by 
the  employed  navigation  transmitters.  These  excited  lower 
hybrid  waves  have  much  shorter  (larger)  wavelengths  (wave 
numbers)  than  those  of  the  injected  VLF  waves.  The 
broadening  introduced  by  this  mechanism  is  characterized 
by  the  fact  that  it  produces  wave  frequency  spectra  that 
exhibit  suppressed  field  intensit y  at  the  carrier  frequency 
and  enhanced  intensity  at  a  discrete  set  of  frequencies  on 
either  side  of  the  carrier  frequency.  Such  spectra  are  less 
frequently  observed  than  those  produced  by  the  nonlinear 
scattering  mechanism.  The  lower  hybrid  waves  excited  by 
this  mechanism  have  enlarged  wave  vectors;  the  apparent 
broadening  due  to  the  Doppler  effect  is  comparable  to  that 
calculated  for  the  first  mechanism. 

We  conclude  that  the  two  suggested  mechanisms  can 
contribute  additively  to  the  observed  spectral  broadening  of 
injected  VLF  waves  reported  by  Bell  et  al.  [1983]  and  others. 
This  spectral  broadening  phenomenon  provides  additional 
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evidence  to  support  the  mechanism  suggested  by  Lee  and 
Kuo  [1984]  for  the  excitation  of  lower  hybrid  waves  by 
injected  VLF  (whistler)  waves. 
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The  combined  operation  of  an  HF  or  MF  ground  transmitter  and  a  VLF 
transmitter  for  enhanced  ionospheric  heating  is  discussed.  The  H  F  or  M  F  transmitter, 
operated  in  a  pulsed  mode,  can  preferentially  produce  short-scale  density  striations 
that  can  render  the  nonlinear  mode  conversion  of  the  subsequently  launched  V  LF 
waves  into  lower  hybrid  waves.  In  addition  to  the  mode  conversion  process,  the 
VLF  waves,  if  intense  enough,  can  also  excite  meter-scale  density  striations  and 
lower  hybrid  waves  via  parametric  instabilities.  Intensified  density  striations  and 
enhanced  airglow  are  expected,  and  they  can  be  detected  by  incoherent  backscatcer 
radars  and  photometers,  respectively.  The  feasibility  and  planning^of  the  proposed 
experiments  are  addressed. 


1.  Introduction 

An  experimental  scheme  is  discussed  for  enhanced  ionospheric  heating  by  the 
combined  operation  of  two  powerful  ground  transmitters.  One  is  an  HF  or  MF 
heater  wave,  and  the  other  is  a  VLF  heater  wave.  The  proposed  scenario  is  as 
follows. 

The  ionosphere  is  illuminated  by  an  HF  or  MF  wave  first  and,  then,  a  VLF 
wave  subsequently  (see  Fig.  1).  It  has  been  known  both  theoretically  and  experi¬ 
mentally  that  short-scale  (i.e..  less  than  a  few  meters)  ionospheric  irregularities  can 
be  excited  within  seconds  by  HF  heater  waves,  while  it  takes  tens  of  seconds  or 
longer  for  large-scale  (say,  hundreds  of  meters  and  longer)  ionospheric  irregularities 
to  be  generated  (see,  e.g.,  Gurevich,  1978;  Fejer,  1979;  Frey  et  al.,  1984;  CosTER 
et  al.,  1985).  As  shown  in  Section  2  of  this  paper,  short-scale  ionospheric  density 
striations  can  effectively  cause  the  nonlinear  scattering  of  VLF  waves  into  lower 
hybrid  waves.  Hence,  we  only  need  to  operate  the  HF  or  MF  heater  in  pulsed  mode 
to  assure  the  excitati  .»n  of  meter-scale  ionospheric  irregularities.  During  the  vertical 
ionospheric  heating,  the  HF  heater  wave  frequencies  should  be  less  than  the  foFi  for 
overdense  heating  of  the  ionospheric  F  region.  When  a  MF  heater  is  used,  the  wave 
frequency  is  required  to  match  the  local  electron  gyrofrequency  in  the  ionosphere  to 
produce  short-scale  ionospheric  irregularities  via  electron  cyclotron  heating  (Lee  ei 
al.,  1986). 
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Fig.  I  The  proposed  experlmem  for  enhanced  ionospheric  heating.  Illumination  of  the  ionosphere  by 
a  HF  or  MF  ground  transmitter  first,  and  then  a  VLF  ground  transmitter  subsequently. 


In  addition  to  the  nonlinear  scattering  process,  the  VLF  waves,  if  intense 
enough,  can  also  excite  lower  hybrid  waves  and  ionospheric  irregularities  via 
parametric  instabilities  (Lee  and  Kuo,  1984).  The  suggested  diagnoses  for  the 
expected  ionospheric  effects  shall  be  discussed  in  Section  3.  In  principal,  the  two 
different  processes  producing  the  lower  hybrid  waves  can  be  experimentally 
distinguished. 
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2.  Heater  Wave-Induced  Ionospheric  Effects 

2. 1  Mode  conversion  via  nonlinear  scattering  process 

Various  plasma  instabilities  can  be  excited  by  HF  heater  waves  during 
overdense  ionospheric  heating.  In  seconds,  short-scale  ionospheric  density  striations 
can  be  generated  concomitantly  with  Langmuir  waves,  upper  hybrid  waves  etc. 
However,  it  requires  tens  of  seconds  for  the  generation  of  large-scale  ionospheric 
irregularities  by  self-focusing  instability  or  thermal  filamentation  instability.  Large- 
scale  irregularities  can  give  rise  to  phase  and  amplitude  scintillation  of  beacon 
satellite  signals.  These  wave  interference  patterns  may  not  be  associated  with 
significant  attenuation  of  the  radio  signals.  By  contrast,  the  short-scale  irregularities 
are  able  to  cause  anomalous  absorption  of  the  radio  signals  via  the  nonlinear 
scattering  process  described  below  that  converts  electromagnetic  wave  energy  into 
electrostatic  wave  energy  (eventually,  plasma  kinetic  energy)  in  the  ionosphere. 

For  simplicity,  a  ducted  whistler  (VLF)  wave  propagating  along  the  i  direction 
is  assumed,  viz.. 


£o  =  (.r  -t-  /v  )£oexp[((Ai);  ~  wi.f)],  ( 1 ) 

Further,  short-scale  ionospheric  irregularities  induced  by  an  HF  or  MF  heater  wave 
are  represented  by 


Sn  =  dnexp(iky).  (2) 

Solving  the  wave  equation,  we  can  show  that  the  ionospheric  irregularities  scatter 
nonlinearly  the  circularly  polarized  (whistler)  wave  into  an  elliptically  polarized 
wave.  In  other  words,  a  linearly  polarized  component  of  the  scattered  wave  is 
introduced  by  the  EM  wave-induced  field-aligned  ionospheric  irregularities. 

The  scattered  whistler  (VLF)  wave  field  can  be  expressed  as 

£,  =  {(.V  -I-  iy)£cp  +  >’£iplexp[/(^oZ  -  wot)],  (3) 

where  £;p=  £b(<5n/  No)  representing  the  wave  field  intensity  of  the  circularly  polarized 
component  while  £ip  represents  that  of  the  linearly  polarized  component,  .Vo  is  the 
unperturbed  plasma  density.  It  is  found  that  £p  and  £ip  are  related  by 


where  A  =  2nl  k  and  Ao=2ntko  are  the  scale  size  of  ionospheric  irregularities  and  the 
wavelength  of  the  whistler  wave,  respectively.  If  a  VLF  wave  at  the  frequency  of  10 
kHz  is  injected,  and  the  ionospheric  conditions  are  assumed  to  be  /c  (electron 
cyclotron  frequency)=  1.3  MHz  and /p  (electron  plasma  frequency)  =  6  MHz,  the 
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whistler  wavelength  in  the  ionosphere  is  found  to  be  about  500  meters.  .According  to 
Eq.  (4),  fip  dominates  over  £cp  provided  that  Hence,  meter-scale  irregularities 

produced  by  HF  of  MF  heater  wave  can  effectively  render  the  nonlinear  mode 
conversion  of  a  VLF  (whistler)  wave  into  a  lower  hybrid  wave.  The  parallel  and 
perpendicular  wave  vectors  of  the  lower  hybrid  wave  are  identical  to  the  wave  vector 
of  the  ducted  whistler  wave  and  that  of  the  field-aligned  ionospheric  irregularities, 
respectively. 

2.2  VLF  wave-excited  parametric  instabilities 

While,  tn  principle,  the  nonlinear  scattering  process  occurs  at  any  field  intensity 
of  the  incident  VLF  wave,  parametric  instabilities  can  be  triggered  by  intense 
enough  VLF  waves  as  suggested  in  Lee  and  Kuo  (1984).  Short-scale  ionospheric 
irregularities  associated  with  zero-frequency  modes  can  be  excited  simultaneously 
with  lower  hybrid  waves  by  the  VLF  waves.  Therefore,  the  HF  or  MF  heater 
wave-induced  short-scale  ionospheric  irregularities  can  be  further  strengthened  by 
the  VLF  heater  waves. 

3.  Diagnoses  of  Expected  Ionospheric  Effects 

It  is  generally  believed  that  spread  F  echos  in  ionograms  and  scintillation 
phenomena  are  caused  by  large-scale  (say,  hundreds  of  meters)  ionospheric  irregu¬ 
larities.  The  generation  of  large-scale  ionospheric  density  irregularities,  therefore, 
can  be  sensed  by  ionosondes  and  scintillation  measurements  of  radio  stars  or  beacon 
satellite  signals.  By  contrast,  the  presence  of  field-aligned  short-scale  (a  few  meters 
and  less)  irregularities  can  be  detected  by  backscatter  radars. 

Backscatter  radar  measurements  can  also  be  made  for  detecting  the  excited 
lower  hybrid  waves  when  the  radar's  beam  angle  is  90°  with  respect  to  the  earth's 
magnetic  field.  The  recorded  incoherent  backscatter  radar  spectrum  should  look  like 
a  double-humped  ion  spectrum  peaking  at  the  lower  hybrid  wave  frequency.  Since 
lower  hybrid  waves  can  accelerate  electrons  effectively  along  the  geomagnetic  field, 
airglow  enhancement  due  to  the  impact  excitation  of  neutrals  by  energetic  electrons 
can  be  expected.  Measurements  of  airglow  at  6300  A  (red),  5577  A  (green),  and  even 
shorter  wavelengths  may  possibly  be  observed.  In  addition,  the  techntque  described 
in  Carlson  et  al.  (1982)  for  finding  the  height  distribution  of  enhanced  plasma  lines 
is  capable  of  distinguishing  further  the  electron  acceleration  by  Langmuir  waves 
and/or  upper  hybrid  waves,  and  lower  hybrid  waves  that  are  produced,  respectively, 
by  HF  and  VLF  heater  waves. 

There  are  several  ionospheric  heating  facilities  operated  at  the  frequencies  of  a 
few  MHz  in  the  HF  or  MF  band  in  the  U.S.A.,  Europe,  and  the  U.S.S.R.  VLF 
transmitters  have  been  used  at  different  locations  for  the  study  of  wave-particle 
interactions  in  the  magnetosphere.  However,  it  may  not  be  possible  to  use  the 
existing  HF  (or  MF)  and  VLF  transmitters  for  the  proposed  experiments.  Since 
VLF  waves  in  a  narrow  frequency  band  can  be  generated  by  lightning  storms,  the 
experiments  can  be  carried  out  with  available  ionospheric  heating  facilities  during 
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lightning  storms.  It  is  not  necessary  to  have  lightning  storms  occurring  nearby.  If 
lightning  occurs  at  the  conjugate  location  in  the  opposite  hemisphere,  the  lightning- 
induced  VLF  waves  can  propagate  in  ducted  whistler  wave  modes  and  bounce  back 
and  forth  for  several  cycles  before  they  die  out. 

This  work  was  supported  by  NASA  contract  NAG5-I055.  Lincoln  Lab  grant  (under  the 
Air  Force  Contract  to  Lincoln  Lab  H  FI9628-80-C-0002)  and  the  RADC  contract  through 
Southeastern  Center  for  Electncal  Engineering  Education  at  the  Massachusetts  Institute  of 
Technology,  and  by  the  NSF  grant  ATM-8713217  and  the  AFOSR  contract  .•kFOSR  88-0127 
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RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
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